Does the Big Rip survive quantization? 
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m t It is known that certain quantum cosmological models present quantum behavior for 

large scale factors. Since quantization can suppress past singularities, it is natural to in- 

O 

£N) ■ quire whether quantum effects can prevent future singularities. To this end, a Friedmann- 

• Robertson- Walker quantum cosmological model dominated by a phantom energy fluid is 

D ■ 

Ph ' investigated. The classical model displays accelerated expansion ending in a Big Rip. The 
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X 



quantization is performed in three different ways, which turn out to lead to the same result, 
namely there is a possibility that quantum gravitational effects could not remove the Big 



(N 

> 1 Rip 

00 . 
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' I. INTRODUCTION 

cr 

si 

bJT), Ever since there appeared compelling evidence that the Universe is expanding at an increasing 

n 

rate many ideas have been put forward to make sense of this unexpected discovery. One of the 



most popular is the assumption that the Universe is homogeneously permeated with a mysterious 
dark energy which, although unseen, dominates the present energy content of the Universe (for 
an up-to-date review see [4]). Accelerated expansion driven by negative pressure is characteristic 
of a wide variety of proposals such as cosmological constant {^quintessence (4], /c-essence Q], 



a, 



braneworlds Chaplygin gas [2], quintom [8j|, and holography [9]. The generalized Chaplygin 



gas, described in terms of a complex scalar field, has been first proposed to unify dark matter 
and dark energy in and Quartessence [lj], which also relies on exotic equations of state, 
and a dark fluid implemented by a complex scalar field with an exponential potential [13|] are 
other attempts at unifyingdark matter and dark energy. A striking candidate to dark energy is 
a so-called phantom field [lj] whose equation of state is p = wp with w < — 1. A recent analysis 
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[151 ] shows that the possibility that our Universe contains phantom energy cannot be ruled out 
by the present observational data. Instead of decreasing as the Universe expands, the phantom 
energy density grows without bound and eventually dominates any other form of energy. The 
scale factor, as well as the phantom energy density, becomes infinite a finite time from now. This 
hypothetical catastrophic end of the Universe, with everything torn apart by the infinite phantom 
energy density, has been christened a "Big Rip" [it] ]. 

Several suggestions have been made either to shun phantom energy or to retain it while avoiding 
the Big Rip. The most radical of such proposals do away with dark energy altogether either by 
ascribing the observed acceleration of the Universe to back-reaction of cosmological perturbations 
\v\ or to a postulated modification of _general relativity through the addition of inverse curvature 
terms to the Einstein-Hilbert action [18]. Other models without dark energy are based on such 
disparate schemes as taking the Universe to be an expanding spherical 3-brane [19y or imputing 



20]. 

2, 



the origin of the cosmic acceleration to the short range interaction of fundamental particles 
Cosmological evolution in the presence of phantom energy without a Big Rip seems possible 
as in the so-called "hessence" models [22j, in which a non-canonical complex scalar field plays the 
role of dark energy and the equation-of-state parameter can cross the phantom divide w = — 1 . It 
has also been claimed that avoidance of the Big Rip in phantom cosmology can be accomplished 
by gravitational back-reaction 23]. 



Phantom-like behavior can arise in a dilatonic brane-world scenario with induced gravity 



or be caused by a non-Hermitian but CPT symmetric Hamiltonian 



251 ]. At the phenomenological 



evel, in which phantom energy is described as a perfect fluid, thermodynamical considerations 



261 ] suggest that phantom energy has negative entropy, which appears to defy a sensible physical 
interpretation. 

There have been studies not only of classical cosmological scenarios motivated by quantum 
effects and containing phantom energy [27] but also of the influence of quantized fields on the 
evolution of classical cosmological models containing a Big Rip [28], in addition to attempts at 
reformulating semiclassical effects in terms of an effective phantom fluid in loop quantum cosmology 



291 ]. There are some examples of quantum cosmological models that present quantum behavior 
for large scale factors jscj. Since, in many cases, quantization has the power to suppress past 
singularities, it is reasonable to inquire whether quantum gravitational effects might prevent the 
Big Rip, which is a future singularity. 

A series of investigations of quantum effects in phantom cosmology based on one-loop effective 
action corrections or back-reaction of quantum fields indicate that the future singularity may be 



3 



softened or even suppressed 



331 ] . It has recently been found that modifications to Friedmann's 



equation induced by loop quantum gravity lead to avoidance of certain future singularities, in- 
cluding the Big Rip [34]. Similarly, a quantum phantom cosmology based on the Wheeler-DeWitt 
equation and the representation of phantom energy by a scalar field has indicated that wave packets 
tend to disperse in the region that corresponds to the Big Rip, thus removing the future singularity 



35]. 



Since there are several independent ways of doing quantum cosmology, and it is not at all clear 
how they are related, one should try to find out what these different formulations have to tell about 
future singularities. One widely explored approach to quantum cosmology consists in describing 
the matter content by means of fluids, making use of some action principle to describe gravity plus 
matter and then performing a canonical quantization. To our knowledge, within this formalism no 
analysis of quantum cosmological models containing phantom energy, with emphasis on the possible 
avoidance of the Big Rip by quantum gravitational effects, has appeared so far. As a preliminary 
step in this direction, here we investigate a quantum cosmological model dominated by phantom 
energy, the latter mimicked by a perfect fluid. As will be seen, this quantization scheme is worth 
pursuing since it is exactly soluble and produces interesting results. Of course it would be more 
realistic not to assume that the phantom fluid dominates the dynamics of the Universe throughout 
its history, as we do to make the model tractable. Nevertheless, this may be considered not too 
serious a weakness since we are interested in the possible influence of quantum effects on the very 
late time dynamics. 

The quantization is performed in three different ways, in order to reach a more satisfactory level 
of generality and gain confidence that the results are not a mere artifact of a specious method of 
quantization. We find that for a certain class of quantum states a phantom fluid brings about a 
Big Rip, just as it does classically. 

Here is a summary of the paper. In Section 2 the classical model is formulated on the basis 
of the canonical formalism of Schutz and the classical equations of motion are solved. In the case 
of a phantom fluid the model displays a Big Rip. In Section 3 a quantization is performed with 
a particular order of noncommuting operators in the Wheeler-DeWitt equation, which requires a 
special inner product. The expectation value of the scale factor grows without bound and becomes 
infinite at a finite value of cosmic time. In Section 4 the order of noncommuting operators is left 
arbitrary so that the quantization can be carried out with the standard inner product. In Section 
5, after reduction of the classical Hamiltonian to that of a free particle with the help of a canonical 
transformation, a new quantization is achieved with the standard inner product. The results of 



4 



Sections 4 and 5 confirm those of Section 3: accelerated expansion and a Big Rip. Section 6 is 
devoted to our conclusions. 

II. THE CLASSICAL MODEL 

n 

Our starting point is the canonical formalism developed by Schutz [3l( to describe gravity plus 
a perfect fluid. The action is 

S=[ A\/^E+2 / d 3 xVhh ab K ab + [ d^x^gp , (1) 
JM JdM JM 

where K ab is the extrinsic curvature and is the metric of the boundary dAi of the four- 
dimensional domain j\A. Units are chosen such that 16ttG = 1, c = 1, and h = 1. The pressure 
p is related to the energy density p by the equation of state p = wp where w is constant. The 
four-velocity of the fluid is expressed in terms of five velocity potentials e, (, (3, 6, S as 

U v = -(e, u + C(3,u + eS, u ) , (2) 
H 

where p is the specific enthalpy, and the comma denotes partial derivative. In the case of 
Friedmann-Robertson- Walker (FRW) models the only nonvanishing potentials are e, 9, and the 
specific entropy S. The normalization condition U U U U = — 1 determines p in terms of the poten- 
tials. 

The FRW metric is 



ds 2 



-N 2 (t)dt 2 + a 2 (t)a ij dx i dx j 



(3) 



where N is the lapse function, a is the scale factor, and aij is the metric of the constant-curvature 
spatial sections. With the help of thermodynamical considerations, and taking the constraints into 
account, insertion of (J3j) into ([!]) yields the reduced action [32] 



{dp a + tpT - NH)dt 



(4) 
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with the super-Hamiltonian 

2 

H = -^--6ka + ^f- . (5) 
24a a 3w y ' 

Here T is a new canonical variable [32] that describes the only remaining degree of freedom of 
the fluid, while k = 0, 1 or —1 depending on whether the curvature of the spatial sections is zero, 
positive or negative. From now on we shall discuss only the flat case (k = 0). 
With k = the equations of motion are 



d , Ar -,v N Pa . d Np 2 a N PT f . 

a = Wa {NH) = " W ' Pa = -d-a {N7i) = "24^ + 3W o^ ' (6) 



supplemented by the super-Hamiltonian constraint 



(> (m = ^L , Pr = -L(NH)=0 , (7) 



2 

24a a 3 ™ v ; 



Solving this constraint for p a we get 



Pa 



where pt is a positive constant, and the negative sign has been selected in order to produce an 
expanding universe. In the cosmic-time gauge N = 1 the first of equations ([6]) becomes 

. = x/2p o _ (3u , + i)/2 ( (10) 

which is solved by 



7 3(t«+l)/2 _ _3(«+l)/2 /PT 



+ A /||3( U ; + l)(t-to) , (11) 



where ao is the present value of the scale factor. If w < — 1 then u; + 1 = — |io + 1| <0 and 
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a(t) 



3(w+l)/2 




3\w + l\(t- t ) 



-2/3|uH-l| 



(12) 



which becomes infinite at the finite future time 



24 a, 



t = to + \ — 



3(tu+l)/2 




Pt S\w + 1\ 



(13) 



This is the famous Big Rip 



III. QUANTIZATION WITH A SPECIAL INNER PRODUCT 

The Wheeler-DeWitt equation results from making the correspondence p a — > p a = —id/da, 
Pt —> Pt = id/dr and requiring = 0, where TC is the operator corresponding to the super- 
Hamiltonian 7i. This means that the time has been chosen as r = — T. With a particular choice 
of operator ordering the Wheeler-DeWitt equation becomes 



i a 



1-31D 



a7 



1 d 2 ^> 



24 da 2 



(14) 



As discussed in 



321 ]. this equation takes the form of a Schrodinger equation id^>/dr = with a 



Hermitian Hamiltonian operator H if the inner product is 



POO 

($,¥)=/ a 1 - 3 ™ <Z>(a,T)*y(a,T)da 
Jo 



(15) 



By superposing stationary states a solution to (|14p can be found 32| in the form (7 is a positive 
constant) 



VP(a, t) = a exp j 



.3(1-10) 



4[ 7 -i3(l-u;) 2 r/32] 



(16) 



The expectation value of the scale factor is 



3J 



(a)(r)«(6 2 + r 2 ) 1/3M) , b 



3(1 -wf 
32 7 



(17) 
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Note, however, that r is not cosmic time. Recall that r = —T, where T is the only remaining 
degree of freedom of the fluid. It follows from the first of equations ([7]) that \N\ = a 3w . Since 
(a) plays the role of the scale factor in the quantum theory, we write dt = (a) 3w dr. Therefore, 
choosing a common origin for both times, 

toe f (b 2 + A 2 r /(1 ~ w) dA . (18) 
Jo 

Near infinity the above integral behaves as 



X 2w/(l-w) dXocT (l+w)/(l-w) _ ( 19 ) 

If w < —1 the integral is convergent and t — > t r i p for r — > oo, where t r i p is finite. Thus the 
expectation value of the scale factor becomes infinite at the finite cosmic time t r i p which corresponds 
to r = oo. The quantum model displays a Big Rip just like the classical model. 

One might suspect that this result lacks generality in view of: (i) the particular order of non- 
commuting operators employed in the construction of the Wheeler-DeWitt equation (|14|) ; (ii) the 
unusual inner product (|15p that, oddly enough, depends on the equation of state through w. We 
proceed to address these issues. 

IV. OPERATOR ORDERING AND QUANTIZATION WITH THE STANDARD INNER 

PRODUCT 



The classical Hamiltonian corresponding to the Wheeler-DeWitt equation (114p is 



,3lD-l 



H = -^pi, (20) 
to which we associate the Hamiltonian operator 

1 1 

H^ v = — -(a x p a a' i p a a u + a u p a a'"p a a x ') , \ + [i + v = '5w-l . (21) 

This operator is formally self-adjoint for any choice of the ordering parameters A, /i, v with the 
standard inner product 



/•OO 

($,*) = / ^(a,T)*^(a,T)da 
Jo 



(22) 
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Now the Wheeler-DeWitt equation takes the form 

= H X ^ , (23) 
whose stationary solutions ^>(a,r) = i\){a)e~ %Er satisfy 

H x ^4> - Ei, = . (24) 
After some algebraic labor this equation can be put in the form 



a 2 ^ + 2ra^- + [e + 24£a 2 ( 1 "% = (25) 
cter da 

where e = [v{2r — A — 1) + A(2r — v — 1)] and r = (3w — l)/2. The solutions of this equation are 



36] 



^a) = a^)/* Z (^l a l-r\ r/1 (26) 



'U-r 



1 - 2r x 2 



where 



p = —\l{—) - < 27 » 

and Z p is a Bessel function either of the first or of the second kind of order p. We shall consider 
only Z p = J p , the Bessel function of the first kind, in which case a complete solution to equation 
([231) is given by 



* E (a, r)=A e~ iE - a^)/2 Jp (^-a 1 ^ 



(28) 



Our next step is to generate a normalizable wave packet by superposition of the solutions f|28[) in 
the form 



Setting A = V24£/(l - r) this takes the more convenient form 



(29) 
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tf(a,r) = /°°dAC(A)e-^ 1 - r ) 2A2T / 24 a( 1 - 2r )j p (Aa 1 - r ) 

JO 



(30) 



The choice 



C(A) = A p+1 e" 7A2 , 7 >0 



allows the integration to be performed in closed form and we obtain [3 



(31) 



^(a,r) = a 



-1/2 



,1-r 



2 7 + iii#i 



i+p 



cxp 



2(l-r) 



47+ 6 ■ 



The expectation value of the scale factor is defined by 



(a)(r) 



J™a\V(a,T)\ 2 da 
/ °°|*(a,r)|2da 



and a straightforward computation furnishes 



(32) 



(33) 



(a)(r) 



2(l-r)(l+p)+l 
2(l-r) 



1 

2^ 



,2 , (l-r) 4 . 



r(p + i) 

Recalling that r = (3io — l)/2, equation (j3~4"|) yields 



7 + 



576 



l/2(l-r) 



(34) 



(a)(r)oc(6 2 +r 2 ) 1 /3(i-) , (35) 

the same as (|17p . Thus, an arbitrary operator order with the standard inner product leads to the 
same result as the one obtained when a special operator order is adopted which requires the strange 
inner product ()15[) . This lends support to our previous claim that there exist evolving quantum 
states such that the phantom energy gives rise to accelerated expansion and a Big Rip. 

It should be stressed that we are trying to find out whether quantum effects necessarily remove 
the future singularity. If this were the case, for any quantum state the expectation value of the 
scale factor would have to remain finite for any value of the cosmic time. We have just presented 
an example of a quantum state such that the expectation value of the scale factor becomes infinite 
in a finite cosmic time. This shows that, at least in this model, and for the above quantization 
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procedure, it is not generally the case that quantum effects prevent the Big Rip. This reasoning 
also justifies our particular choices (I28p and (I3ip . which allow the integration to be performed and 
the wave packet to be expressed in terms of elementary functions, which, in turn, makes easy the 
computation of the expectation value of the scale factor. 



V. CANONICAL TRANSFORMATION AND ANOTHER QUANTIZATION WITH THE 

STANDARD INNER PRODUCT 



The phase-space transformation (a,p a ) — > (x,p) defined by 



2\/l2 ,(l-«,)/2 « ,3 - 1,/2 



-Ki-^ " " = ^7ij-* (36) 

is canonical inasmuch as the Poisson bracket {#,p}( a ,p a ) = 1- In terms of these new canonical 
variables the classical Hamiltonian (I20B reduces to the free-particle Hamiltonian 



# = y, (37) 



and the corresponding Wheeler-DeWitt equation takes the form 



1 d 2 V 

2 ~dx T 



(38) 



This is the Schrodinger wave equation for a free particle on the half-line [0, oo). Self-adjointness 



of the Hamiltonian operator requires 



381 ] that the the domain of H be restricted to those wave 



functions that satisfy the boundary condition 



*(0,r) =£*'(0,r) (39) 

where 5 r ' = d^/dx and £ G [0, oo]. For the sake of simplicity we shall consider only the cases £ = 
and £ = oo, that is, either the wave function or its derivative vanishes at x = 0. The propagator 
for arbitrary £ is known [39J, but it is sufficiently complicated to prevent us from constructing 
manageable wave packets. However, inasmuch as we are searching for examples of quantum states 
that lead to a future singularity, if we are successful our particular choices of £ will still allow us 
to conclude that it is not true that the Big Rip is suppressed in general, that is, for arbitrary £. 
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In short, our special choices for £ are technically convenient and do not impair the logic of our 
reasoning. 

The propagators for the cases £ = and £ = oo are given by [39 ] 



Go(x, x' , r) = G(x, x' , r) — G(x, — x', r) , Goo(x, x' , r) = G(x, x' , r) + G(x, — x' , r) , (40) 
where G(x,x',t) is the usual free-particle propagator, which in our case (m = 1, K = 1) takes the 



form 



40fl 



G(x, x' , r) 



27rir 



1/2 



exp 



(a; — a; 
27 



A2 



(41) 



Let us first consider the initial normalized wave function 



vr / 



-/3x 2 



, P = a + i"f , a > . 



(42) 



Since its derivative vanishes at x = 0, it satisfies the boundary condition (|39p with £ = oo. 
Therefore 



* e (x,r)=/ Goo(x,y,T)V e (y,0)dy 



G(x,y,r)^(y,0)dy , 



(43) 



where we have taken advantage of the fact that x I' e (x, 0) is an even function to extend the integration 
to the whole real line. A simple computation yields 



».(*,r) = ( T ) (l + 2i/3r 



-1/2 



cxp 



1 + 2i/3r 



In accordance with (|36p the expectation value of the scale factor is defined by 



(44) 



(a) e (r) = 

and a simple calculation gives 



'3(1 -w)' 


2/3(1— w) r 




Jo 



3.2/3(1-10) |^ e ( X)T )|2 d:r 



(45) 



(o) e (r) = -= 

V 7T 



3(1 



w 



2V12 



2/3(l-w) 



5 — 3w 



3(1 - w) 



4d 2 r 2 + (1 - 2yr)< 
27 



1/3(1-uj) 



(46) 
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whose asymptotic behavior coincides with that of (|17|) . 
Finally, consider the initial normalized wave function 



V o (x,0) = (^^) 1/4 *e-^' 2 , f3 = a + i 1 , a > , (47) 



which satisfies the boundary condition (|39p with £ = 0. As in the previous case, the fact that 
fy ( x , 0) is an odd function allows us to write 



roc roo 

V {x,t)= G (x,y,T)V o (y,0)dy= G(x,y,T)V o (y,0) dy , (48) 

JO J-oo 

which gives 



^o(x,t) 



/ 128a 3 \ 1/4 



V vr 



7 



1 + 2i/3r)" 3 / 2 x exp 



1 + 2if3r_ 



(49) 



The expectation value of the scale factor is now given by 



(a)o(r) = ^ 
or, equivalently, 



3w 



"3(1 



w 



2/3(l-tu) 



W 



.3(1 



3w 



w 



4a 2 r 2 + (1 - 2 7 r) 
2^ 



2-1 l/3(l-w) 



(50) 



(a)o(T) = ^) <a)e(T) • (51) 

Thus, we recover the same asymptotic behavior of the previous case, which entails accelerated 
expansion and a Big Rip. 



VI. CONCLUSIONS 



We have investigated the quantum features of a FRW cosmological model dominated by a 
phantom energy fluid. The classical model was formulated on the basis of the canonical formalism 
of Schutz, which takes into account the degrees of freedom of the fluid. In the case of a phantom fluid 
(equation of state p = wp with w < — 1) the classical equations of motion predict an accelerated 
expansion for the Universe ending in a catastrophic Big Rip a finite time from now. Then the model 
was quantized in three different ways. Firstly, the Wheeler-DeWitt equation was set up with a 
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special choice of the order of noncommuting operators, which requires an unusual inner product 
that depends on the equation of state for the fluid. It was found that the expectation value of the 
scale factor becomes infinite at a finite cosmic time, and the Universe winds up in a Big Rip. Next, 
an arbitrary order of noncommuting operators was allowed for the purpose of quantizing the model 
with the standard inner product. A particular wave packet was constructed, and it was found that 
the expectation value of the scale factor reproduces the previous behavior. Finally, a third method 
of quantization was explored based on the previous reduction of the classical Hamiltonian to that 
of a free particle by means of a canonical transformation. Since the propagator of the Schrodinger 
equation for a free particle on the half-line is known, it is possible, in principle, to find explicitly the 
wave function at any time once its form is given at an initial time. It turns out that two initial wave 
functions obeying different boundary conditions give rise to essentially the same expectation value 
of the scale factor with exactly the same behavior predicted by the two other quantization schemes. 
This indicates that our results are not merely an artifact of a spurious quantization procedure. The 
quantum behavior of the scale factor appears to be robust, and supports our conclusion that, at 
least in our model, quantum effects do not, in general, prevent a Big Rip. 

We do not claim that our result is conclusive. In order to check its generality, the matter 
content should be enriched to describe the early universe more realistically, which we intend to 
do in a future work. Since the phantom energy dominates the late time dynamics, we suspect 
that the enrichment of the matter content will not change our result, but this conjecture must be 
checked. We anticipate that this will not be an easy task because it will surely be much more 
difficult to obtain exact normalizable solutions to the Wheeler-DeWitt equation. We do not expect 
that approximate solutions, such as those furnished by the WKB method, are suitable to discuss 
the existence of quantum singularities, that is, we believe that only a fully quantum treatment can 
give us a solid clue as to the persistence of future singularities in the quantum regime of phantom 
cosmological models. We also feel that a study of the Bohmian trajectories of the scale factor is 
likely to be helpful in clarifying whether quantum gravitational effects can suppress the Big Rip. 
This is an investigation we expect to undertake in the near future. 

Our result differs from those obtained by means of other formalisms 



33, 



34, 



351 ] . namely that 



the Big Rip is softened or avoided by quantum effects. On the one hand, since it is not known 
how these different formulations of quantum cosmology are related, we confess to our ignorance 
about the reason for the discrepancy. On the other hand, it is hard to make clear-cut statements 
on presence or absence of cosmological singularities in the quantum regime because there is no 
general agreement on what constitutes a quantum singularity. These issues certainly deserve a 
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deeper investigation. 
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